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An infrared study of the phonon spectra of ZrW^Os as a function of temperature which includes 
the low energy (2-10 meV) region relevant to negative thermal expansion is reported and discussed 
in the context of specific heat and neutron density of states results. The prevalence of infrared active 
phonons at low energy and their observed temperature dependence are highly unusual and indicative 
of exotic low-energy lattice dynamics. Eigenvector calculations indicate a mixing of librational and 
translational motion within each low-frequency IR mode. The role of the underconstrained structure 
in establishing the nature of these modes, and the relationship between the IR spectra and the large 
negative thermal expansion in ZrW^Os are discussed. 



The tendency of solids to expand when heated is one of 
the most pervasive and familiar phenomena in solid state 
physics 0, however, there are some compounds in which 
the opposite phenomenon occurs 0,0. In ZrM^Os this 
negative thermal expansion (NTE) is large, isotropic and 
persists over a broad temperature ranged IE IE IS IE E| • 
These characteristics enhance our interest in both the 
fundamental physics of ZrW^Os, and the potential ap- 
plications of such a substantial NTE effect, which are 
widespread Q. A defining feature of ZrVV^Os is the ex- 
istence of a 1-fold coordinated oxygen site for each WO4 
tetrahedron. This terminal, or unconstrained 0, oxygen 
creates a structural openness along the high-symmetry 
(111) axes and is thought to be influential in the low- 
energy dynamics and crucial to the phenomenon of NTE 
in ZrW 2 8 . 

Thermal expansion phenomena generally come from 
anharmonic phonon dynamics 0, IE IE 0- l3 For 
ZrW208 analysis of temperature denendence|7 Isl Ifll Hoj] 
indicates that the low energy region should be of primary 
importance. To probe the low-energy phonon dynamics 
we use temperature dependent infrared spectroscopy cov- 
ering the range from 5000 cm -1 to 16cm" 1 (600-2meV), 
which encompasses all the optic phonons of ZrW20§. We 
find that the phonon-related peaks in the infrared spec- 
tra extend to unusually low energy (~ 3meV), and that 
there is strong and unusual temperature dependence in 
the low-energy region, reflecting evidence of unconven- 
tional and anharmonic behavior. We suggest that under- 
standing the mechanism of NTE and understanding the 
exotic nature of the infrared spectra from ZrV^Os are 
two essentially similar problems, both rooted in exotic 
low energy dynamics. 

One of the central unresolved issues for ZrW20g has 
to do with the nature of the eigenmode(s) responsible 



for NTE. While earlier work has emphasized the role of 
transverse oxygen vibration0JH IE m El f° r the mech- 
anism of NTE, Cao et al.|ll|,ll2j have proposed a model, 
based on XAFS data, in which translation of WO4 tetra- 
hedra along (111) directions plays an important role. Us- 
ing eigenvector calculations together with our data, we 
infer that transverse oxygen motion, in the form of libra- 
tion, tends to be accompanied by translation within each 
low-energy mode (thus they are essentially inseparable). 
This mixing arises due to the unconstrained oxygen and 
the associated openness of the structure along the high- 
symmmetry (111) direction axes, and it appears to be 
essential to the nature of the low-frequency IR spectra. 
This mixing may also be crucial to the mechanism of the 
NTE: while the transverse oxygen motion provides the 
thermal motion that contracts the lattice, the presence 
of a translational component in each low-energy mode 
provides a source of frustration that inhibits lattice in- 
stabilities that would destroy NTE. 

The single-crystal samples are grown using a layered 
self- flux technique^!. Our measurements of the spec- 
ular reflectivity cover the frequency range from 15 to 
5000 cm -1 and the temperature range 20 to 300 K. The 
optical conductivity, (t(u>), is determined from the reflec- 
tivity via a Kramers-Kronig transform with carefully cho- 
sen upper and lower terminations. For the acentric P2\3 
structure of ZrW20s one expects a total of 132 modes: 
3 acoustic modes, 32 triplets of IR optical phonons, 11 
doublets and 11 singlets, associated with the irreducible 
representations 33T+ 11E+ llA ^^Mjjk. O ur results 
are broadly consistent with earlier work |l4Hl5[ll^ . lrRll^ | 
for the higher frequency parts of the phonon spectrum, 
and complementary in our emphasis on the low frequency 
range. 

Figure la shows the real part of the optical conduc- 



2 




200 w ( cm -i) 600 800 1000 




20 30 40 80 85 90 95 100 
CO (cm" 1 ) CO (cm" 1 ) 



FIG. 1: (color online) a) The real part of the conductivity, 
ai(ui), is shown for the entire optic phonon region. The in- 
set shows the low frequency region for temperatures of 20, 
80(dotted) and 240 K. b) and c) show cti(lj) for T=20, 60, 
80, 100, 120, 140, 160, 180, 200, 220, 240, 260 and 300 K 
in two low u> regions, along with neutron density of states in 
(b)0 



tivity, <7i(u>), vs frequency. The infrared active optical 
phonon modes, which generally appear as peaks in a\{w), 
are sufficiently dense that they tend to merge together 
into clusters in ZrW20s- The upper cluster is associated 
with bond stretching motion (including the modes of 
the WO4 tetrahedra) and the middle cluster (180-400 
cm -1 ) with bond bending modes (e.g., ^2 modes). Below 
these two clusters there are additional infrared features, 
shown in the inset, which we associate with librational 
and translational motion of WO4 tetrahedra, as discussed 
below. Unlike the generic modes of the upper clusters, 
the modes in this region are specific to the ZrW^Os struc- 
ture, and it is in this lower region that one expects to find 
modes relevant to the mechanism of NTE. 

As shown in Fig. lb) and c), the peaks in this re- 
gion extend to unusually low energy and exhibit strik- 
ing temperature dependence. They tend to sharpen and 
grow in strength as the sample is cooled; they exhibit 
swept, asymmetric lineshapes, and in some cases substan- 



tial dependence of peak frequency on temperature. The 
largest peak frequency shifts are observed for the broad 
28 cm -1 mode, which shifts from about 27 to 29.5 cm -1 
(+8%), and the 88 cm -1 mode, which shifts from about 
91.5 to 86cm" 1 (-6%). These shifts are unusually large 
and the substantial accretion of spectral weight (presum- 
ably coming from the higher frequency) is beyond the 
scope of conventional phonon models. The pervasiveness 
of strong temperature dependence throughout this low 
frequency region reflects the exotic character of the low 
energy dynamics which is clearly evident in the IR spec- 
tra. 

It is informative to consider the relationship between 
our IR data and specific heatQ. For the fit shown in Fig- 
ure 2a, a lowest energy Einstein mode at 28.5 cm -1 with 
a spectral weight of 2.6 oscillators per unit cell provides 
a good fit to the leading edge of C(T)/T 3 0. This fit 
is highly constraining with regard to the energy of the 
lowest Einstein mode; attempts to fit C{T) with a lower 
optical mode generate excess specific heat at low tem- 
perature as shown in Figure 2a. We therefore conclude 
that the 20 cm -1 peak does not correspond to the zone- 
center endpoint of an ordinary optical phonon branch 
and that the broad 28 cm -1 peak is the signature of the 
lowest optical phonon of ZrW20s- The total spectral 
weight in the peaks between 28 and 100 cm" 1 is in agree- 
ment with the results of our phonon eigenvector calcula- 
tions, discussed below, and the peaks in this region are 
in reasonable correspondence with zone-center intercepts 
of optical phonon dispersion relations calculated for the 
isostructural compound HfWoQgflfl l2fT| as well as neu- 
tron DOS peaks 8J. The 20cm" 1 peak differs from the 
others in that it diminishes substantially below 80 K (the 
others grow) and is not reconcileable with neutron DOS 
and C(T). These two features may suggest a low-density 
defect with an extremely large dipole matrix element. 
Further theory is needed, in conjunction with controlled 
defect studies, to ascertain the origin of this marked ab- 
sorption, and its possible connection with the unusual 
lattice dynamics of the bulk compound. 

Focusing on energetics, one can let the infrared spectra 
provide a context in which to explore the origins of the 
temperature dependence of the lattice parameter, a(T). 
Similar to earlier work0 0, we approach this by asso- 
ciating negative Gruneisen parameters (7) with Einstein 
mode frequencies. We calculate a(T) using a bulk modu- 
lus of 72 GPa [il| and assuming triplet weight 0. Start- 
ing with the frequency of the lowest optic mode from our 
data (at 28cm _1 =3.5mey), we find that associating a 
large negative Gruneisen parameter to this mode alone 
does not allow a good fit to a(T), however, a good fit can 
be obtained simply by including a second Einstein mode 
at higher frequency, as shown in Figure 2b) for 7=-20 at 
28 cm -1 and 7=- 7 at 88 cm" 1 . The lower energy corre- 
sponds to the lowest energy optical phonon in ZrW^Og, 
which is so low it cuts through the acoustic phonon dis- 



3 




FIG. 2: (color online) a) A fit to the specific heat with a lowest 
infrared mode at 28.5 cm -1 (3.5 meV) is shown by the solid lie 
(through the specific heat data). The dashed curves show the 
excess specific heat generated by failed attempts to fit C(T) 
with a lowest optic mode at 20 or 24 cm -1 , respectively, b) 
a(T) data (open squares) is shown along with a fit based on 
a two Einstein mode model. 

persion and may thereby acquire mixed character [2^; the 
higher energy corresponds to the strongly temperature 
dependent peak near 88 cm -1 (Fig lc). While this 2- 
mode approach provides the minimal context in which a 
good fit can be obtained, one can also do more nuanced 
fits involving additional modes. One can use Griineisen 
parameters for modes above 40 cm -1 from pressure de- 
pendent Raman measurements of Ravindran et a.1. |l5| 
and add an additional term for the 28 cm -1 mode. In 
that more nuanced approach good fits to a(T) can be ob- 
tained when the 7's of the two lowest optic modes add 
to about -20. 

Mechanisms of NTE based on transverse thermal mo- 
tion of oxygens in two-fold coordination (W-O-Zr) have 
been widely considered^ El EJ S El however, Cao 
et al. 0, 0| have proposed an alternative model, based 
on XAFS data, in which translation of WO4 tetrahcdra 
along the high-symmetry (111) axes plays a key role. 
Ultimately the discussion of mechanism centers on the 
eigenvectors of the relevant phonon modes. Toward this 
end, we have calculated the eigenvector for each k=0 
optical phonon mode using a mass-spring model of a 44 
atom unit cell with nearest-neighbor stretching and bend- 
ing interactions and periodic boundary conditions. With 
appropriate bond strengths, the eigenvectors exhibiting 
bond stretching and bending motion are associated with 
high (700-900 cm" 1 ) and moderate (150-400 cm' 1 ) fre- 
quency eigenvalues respectively, as expected. For the 
range below about 120 cm -1 our calculation generates 27 
additional optic modes; of which 21 are associated with 
triplets. The two lowest energy modes are triplets at 32 



and 43 cm -1 respectively, in reasonable correspondence 
with the data. The eigenvectors for these modes exhibit 
a mixture of librational and translational motion, as il- 
lustrated in Figure 3a. A suitable choice of eigenvector 
within the triplet degeracy manifold is required to clarify 
the motion in this way. For the mode shown, clock-wise 
rotation about the (111) axis is accompanied by down- 
ward translation of the two tetrahedra along the (111) 
axis [2^ 

This (111) axis pair is connected, via Zr-centered oc- 
tahcdra, to other tetrahedra as illustrated in Figure 3b. 
The are a total of eight WO4 tetrahedra per unit cell 
arranged in pairs along (111) axes and offset from each 
other in a "spiral staircase" pattern as shown. Motion 
on the other (111) is not the same|27) and the motions 
are correlated in complex ways. 

This mixing of translational and rotational (libra- 
tional) motion is intimately related to the undercon- 
strained nature of the ZrW20s structure, in which each 
tertrahedron has an unconnected (terminal) oxygen at 
its (111) axis apex0- Because there are no strong sec- 
ond bonds for this unconstrained oxygen, translations of 
WO4 tetrahedra along (111) axes do not involve signifi- 
cant bond compression and therefore project to very low 
frequency. There they tend to mix with the low frequency 
librational motion, and it is through this mixing that 
these low energy modes can acquire a dynamic dipole mo- 
ment (since the WO4 tetrahedra carry net charge). Our 
observation of non-zero dipole moment for many low en- 
ergy modes indicates that this mixing is pervasive. Thus 
the richness of the low frequency infrared spectra is con- 
nected to the existence of the unconstained oxygen and 
the related openness along the high symmetry (111) di- 
rections, which is a defining feature of the ZrW20s struc- 
ture and critical to NTE@,|l|g|,[l(J|22 . 

Regarding the origins of anharmonicity in ZrW20s, we 
note that, in the context of simple harmonic picture, the 
uncertainties in the circumfrential oxygen position would 
be substantial. Using the frequency of the lowest optic 
mode (28 cm" 1 ) one gets a theoretical uncertainty in the 
ground state of A9 = y/2h/m r 2 uj ~14° with a cor- 
responding circumferential uncertainty A(r0) ~ 0.4 A. 
These increase to about 40° and 1.2 A by room temper- 
ature. Such large values suggest a possible origin of an- 
harmonicity and motivate consideration of the relevance 
of configurational tunneling or rotation to the low en- 
ergy dynamics. Tunneling, rotation and configurational 
exchange phenomenon have been observed in systems 
with a similar 3-fold rotation axis involving methyl-group 
tetrahedraJU 0] , where the breakdown of the harmonic 
approximation for low-energy librational motion has been 
carefully studied. For methyl-group tetrahedra such as 
CH3I the lowest librational mode frequency is typically 
80 cm -1 , about 3 times higher than Z1AV2O8. That most 
of this difference can be accounted for by the hydrogen- 
oxygen mass difference suggests a roughly comparable 
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FIG. 3: (color online) a) The nature of a low frequency triplet 
mode is shown in terms of WO4 tetrahedra and Zr-centered 
octahedra via a time sequence. As the tetrahedra rotate clock- 
wise, they translate downward along the (111) axis, which is 
shown by the dark line. The two lowest-frequency optic mode 
are similar, but in the other one the (111) axis octahedron 
(white) rotates opposite the tetrahedra. The unconstrained 
oxygen lies on the (111) axis, b) The pair of tetrahedra shown 
in (a) couples to other tetrahedra as shown. There are eight 
WO4 tetrahedra in a unit cell organized in pairs on (111) axes 
(dark lines) which are offset in a spiral-staircase configuration 
as shown. (Illustration by Alison Kendall) 



confining potential to lowest order. A possible differ- 
ence between XCH3 and WO4 tetrahedra arises from the 
nuclear spin, which is ^ for H, and for the O 16 nu- 
cleus, therefore a 120° reorientation of an WO4 tetrahe- 
dron is equivalent to a double permutation among iden- 
tical particles and thus does not produce a distinct state. 
The potential relevance of finite barrier phenomena and 
large oxygen excursions to the the low-energy dynamics 
of ZrW2C>8 remains a question for future work. 

We have provided evidence for highly unusual low- 
energy phonon dynamics in ZrW20s as reflected in our 
infrared spectra. We infer that the lowest optic modes 
in ZrW20g tend to have a mixed librational and transla- 



tional character in which the unconstrained oxygen plays 
an essential role, and it is likely that these modes play a 
central role in the mechanism of negative thermal expan- 
sion. While further work is needed to fully elucidate the 
mechanism of NTE in ZrW20s, the present works shows 
that both librational and translational motions are op- 
erative: the transverse oxygen motion associated with 
libration provides NTE, while translational component 
within each mode frustrates a displacive transition that 
would otherwise remove the unique structural environ- 
ment in which these mixed phonon modes exist. 
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